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HIGHLIGHTS 


•  The  significance  of  reversible  heat  (RH)  is  analyzed  in  a  Li-ion  battery  model. 

•  RH  is  important  in  thin  electrode  cell  and  in  small  solid  particle  cell. 

•  Thin  electrode  cell  and  fine  solid  particle  cell  exhibit  good  thermal  performance. 

•  A  correction  diagram  is  proposed  to  modify  the  results  of  RH  exclusive  models. 

•  RH  is  non-ignorable  for  battery  discharged  at  10  C  or  lower  rates. 
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As  a  typical  heat  source  in  the  operation  of  lithium-ion  (Li-ion)  battery,  reversible  heat  (RH)  makes  a 
considerable  contribution  to  the  temperature  variation  of  the  battery,  and  the  determination  of  this 
contribution  is  of  great  significance,  especially  for  mission-critical  applications.  In  previous  research 
works,  the  determination  of  including  or  excluding  the  RH  was  typically  made  based  on  the  rate  of 
discharge.  This  study,  however,  finds  that  the  importance  of  RH  could  be  more  meaningfully  judged 
based  on  its  contribution  to  the  overall  thermal  response  of  Li-ion  battery  at  various  conditions.  An 
improved  electrochemical— thermal  coupling  model  is  developed  to  incorporate  the  RH.  This  model  is 
applied  to  simulate  the  operation  of  batteries  with  different  physical  properties  (e.g.,  electrode  thickness 
and  active  material  particle  size)  under  different  discharge  rates.  The  simulation  results  are  compared  in 
terms  of  RH  generation  rates,  battery  temperature,  as  well  as  several  other  associated  characteristics.  It  is 
found  that  the  significance  of  RH  to  the  overall  heat  generation  varies  when  different  intrinsic  and 
extrinsic  parameters  are  considered.  This  significance  becomes  more  appreciable  for  batteries  with 
thinner  electrodes  and/or  finer  active  material  particles.  A  correction  diagram  is  then  proposed  to  modify 
the  results  for  thermal  analyses  without  including  RH. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium-ion  (Li-ion)  battery  system  appears  to  be  an  effective 
approach  for  energy  storage  and  is  being  applied  to  address  many 
environmental  and  economic  challenges  in  modern  society. 
Although  numerous  innovations  have  been  made  in  developing 
new  battery  materials,  such  as  high  capacity  loading  anode  [1  ],  auto 
shutdown  coating  and  electrolyte  [2,3]  and  high  voltage  cathode 
[4],  Li-ion  technologies  have  not  yet  been  fully  demonstrated  to 
meet  the  stringent  demand  from  a  variety  of  industries  for  high- 
performance  batteries.  The  reliability  and  safety  issues  arising 
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from  excessive  heat  generation  is  being  considered  as  one  of  the 
most  critical  technical  challenges  inherent  to  Li-ion  battery,  which 
may  lead  to  significant  battery  performance  degradation  or  even 
thermal  runaway  [5  ,  thereby  hindering  the  Li-ion  technologies  to 
achieve  broader  market  penetration. 

In  recent  years,  many  coupling  models  have  been  suggested  to 
solve  the  electrochemical  and  thermal  issues  associated  with  Li-ion 
battery.  For  example,  Gu  et  al.  6]  coupled  an  electrochemical  bat¬ 
tery  model  with  a  thermal  model.  The  electrochemical  character¬ 
istics  (in  terms  of  voltage  and  current  density)  as  well  as  the 
thermal  characteristics  (in  terms  of  heat  generation  rate  and  tem¬ 
perature  variation  under  different  cooling  conditions)  were 
compared  between  both  the  coupled  and  decoupled  models,  and  it 
was  shown  that  an  improved  accuracy  can  be  achieved  when  the 
thermal  model  is  included.  In  Ref.  [7],  the  heat  generation  and 
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Fig.  1.  The  schematic  illustrations  of  a)  the  multi-particle  electrochemical  reaction  cell  model;  and  b)  the  3D  thermal  model  of  the  battery  used  for  validation. 


current  redistribution  in  Li-ion  battery  pack  were  investigated 
based  on  a  coupled  electrochemical-thermal  model  considering 
the  contact  resistance,  and  the  heat  surge  from  inter-cell  connec¬ 
tive  resistance  was  discussed.  However  the  effect  of  the  reversible 
heat  (RH)  in  both  models  was  ignored.  RH  is  also  known  as  the 
entropy  heat  and  is  linked  with  the  entropy  changes  in  the  elec¬ 
trodes  during  the  de-intercalation  and  intercalation  processes.  In 
this  study,  we  call  the  electrochemical-thermal  coupling  model 
that  does  not  include  RH  a  “without  reversible  heat  model”,  or 
“wor-model”.  Whether  the  wor-model  can  accurately  describe  the 
heat  generation  of  a  Li-ion  battery  remains  as  a  challenging  tech¬ 
nical  problem.  Therefore,  a  study  was  conducted  by  Srinivasan  and 
Wang  [8]  on  the  electrochemical  and  thermal  behaviors  of  Li-ion 
battery  using  two-dimensional  unit  layer,  which  indicated  that 
the  RH  is  important  at  all  discharge  rates.  The  experimental  results 
were  also  provided  in  the  literatures  [9-12]  for  demonstrating  the 
entropy  changes  in  different  active  materials.  Jeon  et  al.  [13] 
showed  the  effect  of  RH  to  the  overall  heat  generation  in  Li-ion 
battery  using  LiCo02  and  LiNiCoMn02  as  cathodes  at  different 
discharge  rates  when  the  state  of  charge  equals  to  0.1,  and  it  was 
found  that  the  thermal  response  of  the  battery  with  LiCo02  cathode 
has  a  strong  dependence  on  the  RH.  We  shall  call  the  electro¬ 
chemical-thermal  model  with  the  RH  included  a  “with  reversible 
heat  model”,  or  “wr-model”.  In  order  to  obtain  a  clearer  under¬ 
standing  of  the  heat  generation  across  the  battery,  the  wr-model 
will  be  used  in  this  work.  While  the  previous  research  works 
mainly  demonstrated  the  importance  of  the  RH  on  the  entire  bat¬ 
tery,  this  study  would  consider  the  heat  more  thoroughly  on  a 
component  level,  with  emphasis  on  the  anode  and  the  cathode  in 
an  individual  reaction  cell.  Also,  unlike  the  existing  literatures  that 
only  show  the  impact  of  discharge  currents  on  the  RH  generation, 
this  work  takes  into  consideration  many  other  important  electro¬ 
chemical  factors  such  as  the  electrode  thickness  and  the  active 
material  particle  size  that  may  directly  or  indirectly  affect  the  RH 
generation. 

In  this  study,  the  electrochemical  and  thermal  behaviors  of 
LiMn204-based  battery  are  studied  using  an  improved  coupling 
model  with  the  RH  term  included.  By  changing  the  discharge  rate, 
electrode  thickness  and  particle  size  of  active  material,  the  elec¬ 
trochemical  and  thermal  performances  of  not  only  the  entire 


battery  but  also  the  individual  cell  components  are  investigated. 
For  comparison,  the  simulation  results  based  on  a  wor-model  are 
also  presented.  A  comprehensive  understanding  of  the  significance 
of  the  RH  in  Li-ion  battery  under  various  conditions  is  achieved 
through  in-depth  analysis,  and  a  correction  diagram  that  can 
modify  the  results  for  wor-model  is  thus  proposed,  which  is 
essential  for  future  research  in  battery  heat  generation  and  battery 
thermal  management  system  development. 

2.  Electrochemical-thermal  battery  model 

The  fundamental  governing  equations  of  the  model  in  this  study 
are  mainly  based  on  the  model  developed  by  Smith  and  Wang  [14] 
with  some  modifications  made  by  including  the  double  layer 
capacitance  term  [15]  and  entropy  heat  term.  The  coupling  model 
consists  of  a  ID  electrochemical  part  and  a  3D  thermal  part,  in 
which  the  heat  generated  from  the  electrochemical  part  and  the 
temperatures  obtained  from  the  thermal  part  are  mutually  used  as 
inputs  for  each  other.  The  schematic  of  the  electrochemical  part 
(i.e.,  a  single  reaction  cell  in  lithium  ion  battery)  is  shown  in  Fig.  la, 
and  all  the  calculations  are  performed  in  the  x  direction. 

The  ID  electrochemical  model  has  five  components  from  left  to 
right:  copper  current  collector,  anode,  separator,  cathode  and 
aluminum  current  collector.  During  discharge,  the  lithium  de- 
intercalates  from  the  anode  particles,  resulting  in  an  equal  num¬ 
ber  of  Li-ions  and  electrons  at  the  interface  of  anode  and  electro¬ 
lyte.  The  Li-ions  diffuse  through  the  electrolyte  solution  and 
intercalate  into  the  cathode  active  materials,  while  the  electrons 
flow  through  the  external  circuit  from  the  anode  to  the  cathode  to 
take  part  in  the  electrochemical  reaction.  To  mathematically 
explain  the  electrochemical  reaction  and  the  associated  thermal 
response,  equations  that  describe  the  mass  and  charge  balances  of 
Li-ions  in  the  electrochemical  model  as  well  as  the  energy  con¬ 
servation  in  the  thermal  model  should  be  properly  considered. 

2.1.  t-D  electrochemical  model 

In  the  solution  phase,  the  mass  balance  of  Li-ions  in  the  elec¬ 
trolyte  is: 
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where  the  last  term  on  the  right  hand  side  is  the  contribution  from 
the  double  layers,  note  that  dq_/dq  =  0  and  dq+/dq  =  1  for  a  cation 
adsorbing  electrode  and  dq_/dq  =  1  and  dq+/dq  =  0  for  an  anion 
adsorbing  electrode  [15]. 

The  solid  phase  diffusion  of  the  lithium  inserted  into  the  active 
material  particle  is  governed  by  Fields  second  law  in  the  polar 
coordinate: 


dcs Ds  d  (  20Cs\ 

at  r2  dr  \  dr ) 


The  current  densities  in  the  solid  particle  (is)  and  in  the  solution 
phase  (ie)  are  given  by  Ohm’s  law: 


is 


a  Qs 
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The  relation  between  the  divergence  of  the  current  flow  in  the 
solution  phase  and  the  reaction  current  density  is: 


Vie  = 


The  intercalation  kinetics  in  both  electrodes  can  be  written  in 
the  form  of  the  Bulter-Volmer  equation: 


askc°  5  (Cs5max  —  Cs)  5 


exp 


-  exp 


0.5 F  \ 

~kTv) 


wherein  the  over-potential  q  is  given  by: 


y  =  <Ps-  <pe~u  (7) 

where  the  open  circuit  potential  U  of  the  battery  cell  is  given  by 
Taylor's  first  order  expansion  at  a  reference  temperature  of  25  °C: 


U  =  Ulef  +(T-  Tref)  (8) 

and  the  reference  potential  Oref  of  the  anode  and  the  cathode  at 
25  °C  in  the  equation  is  expressed  elsewhere  [16  .  The  entropy 
change  term  (dLf/dT)  for  cathode  is  expressed  as  [17]: 


=  4.3128  exp(0.5715y)  +  1.2816  sin(-4.9917y) 

-  0.0904  sin(-20.967y  +  12.5788) 

-  0.03135  sin(31.7663y  -  22.4296)  -  4.1453  +  8.1471y 

2 

-  26.0646y2  +  12.766y3  -  0.1843  exp(-^2±g|?) 

(9) 


and  for  the  anode  as  18]: 

du  — 3^1-1  135  exp(-32.963x  + 8.317) 

dT  '  J^1  +749.076  exp(-34.791x  + 8.887)  (10) 

+  0.2698  -  0.852x  +  0.362x2 

where  y  and  x  in  the  Eqs.  (9)  and  (10)  denote  the  ratios  of  the  Li-ion 
concentrations  inside  the  active  materials  during  cycling  and  the 
theoretical  maximal  concentrations  in  the  cathode  and  the  anode, 
respectively.  The  corresponding  fits  are  plotted  in  Fig.  2. 

2.2.  3-D  thermal  model 

Conservation  of  energy  for  a  Li-ion  battery  considering  con¬ 
vection  term  and  heat  generation  term  can  be  expressed  as: 

—  V  •  (fcVT)  +  qr  +  qrev  +  qj  (1 1 ) 

where  the  three  “q”  s  are  the  volumetric  heat  generation  rate  of 
different  heat  sources  from  the  electrochemical  model,  i.e.,  reaction 
heat  (qr),  reversible  heat  (qrev)  and  ohmic  heat  (qj): 


Fig.  2.  Fitting  curves  of  the  entropy  term  (dt//dT)  for  a)  LixC6  from  Eq.  (10);  and  b) 
LiyMn204  from  Eq.  (9). 
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where  l  is  the  thickness  of  the  five  components  (anode,  separator, 
cathode  and  two  current  collectors)  in  the  cell.  In  Li-ion  battery, 
each  current  collector  is  generally  shared  by  two  adjacent  reaction 
cells,  therefore  the  thickness  of  the  current  collector  used  in 
simulation  is  half  of  its  actual  thickness. 

At  the  surface  boundary  of  the  battery,  convective  cooling  is 
considered: 


All  the  parameters  used  in  the  electrochemical  and  thermal 
models  are  listed  in  Tables  1  and  2,  respectively.  The  lv  and  l  in 

Table  2  are  the  thicknesses  of  the  components  in  the  battery  used 
for  validation  and  in  the  battery  model  based  on  Doyle's  work  16], 
respectively. 

3.  Results  and  discussion 


-/<—  =  h(T-T3mb)  (15) 

where  h  refers  to  the  convective  coefficient  and  its  value  is  set  as 
low  as  0.1  W  nrT* 2  k”1.  This  simulates  a  condition  with  a  very  poor 
heat  transfer  rate  (e.g.,  aerospace  applications)  and  can  lead  to  a 
dramatic  temperature  rise  in  batteries,  which  facilitates  a  clearer 
inspection  on  the  differences  in  temperature  elevation  between 
batteries. 

The  density,  heat  capacity  and  thickness  of  each  component  are 
used  to  determine  the  average  values  of  the  density  and  heat  ca¬ 
pacity  of  the  entire  battery: 


pCp  -  ^2 
i 


(16) 


where  k  is  the  thickness  of  the  component  i,  and  L\  denotes  the 
thickness  of  the  entire  battery.  However,  the  thermal  conductivity 
of  the  battery  should  be  considered  separately  in  different  di¬ 
rections.  As  suggested  by  Fig.  lb,  there  exist  series  thermal  resistors 
in  x  direction  and  parallel  thermal  resistors  in  y  and  z  directions, 
therefore  we  define: 


Table  1 

Parameters  used  in  the  ID  electrochemical  model  [16]. 


Parameter 

Negative  electrode 

Separator 

Positive  electrode 

Rs  (pm) 

12.5 

— 

8.5 

Ds  (m2  s"1) 

3.9  x  10  14 

— 

1  x  10"13 

£e 

0.357 

0.724 

0.444 

£p 

0.146 

0.276 

0.186 

£f 

0.026 

0 

0.073 

Cs.max  (m°l  m~3) 

26,390 

— 

22,860 

Cs, 0  (mol  m~3) 

14,870 

— 

3900 

ce  (mol  m-3) 

2000 

2000 

2000 

a  (S  ur1) 

100 

— 

3.8 

0.363 

0.363 

0.363 

Cdi  (F  m-2) 

0.2 

— 

0.2 

Typically,  heat  generation  in  Li-ion  battery  consists  of  three 
sources:  reversible  heat,  reaction  heat  and  ohmic  heat.  Reversible 
heat  (RH),  as  stated  earlier,  accounts  for  the  entropy  changes  in  the 
structure  of  the  electrodes  during  electrochemical  reactions.  Re¬ 
action  heat  refers  to  the  heat  rises  from  the  reaction  current  and 
over-potentials.  And  ohmic  heat  is  mainly  associated  with  the  heat 
generated  during  the  motion  of  lithium/lithium-ions  in  the  solid 
and  solution  phases.  The  generation  rate  of  each  heat  source  is 
related  to  their  respective  amount  of  heat  generated  in  a  certain 
time  range,  and  the  contribution  of  each  of  them  to  the  overall  heat 
generation  gives  an  indicator  of  the  importance  of  the  heat.  With 
regard  to  the  RH,  its  generation  rate  can  be  related  directly  to  the 
discharge  rate,  and  its  contribution  is  relied  heavily  on  the  amounts 
of  ohmic  heat  and  reaction  heat  generated,  which  are  found  to  be 


Table  2 

Parameters  used  in  the  3D  thermal  model  [19—21]. 


Parameter 

Cu  foil 

Negative 

electrode 

Separator 

Positive 

electrode 

Al  foil 

lv  (pm) 

15 

37.5 

35 

24 

20 

l  (pm) 

14 

100 

52 

174 

20 

k  (W  itT1  K"1) 

398 

5 

0.334 

5 

238 

P  (kg  m-3) 

8930 

2500 

1009 

1500 

2702 

C„  o  kg-1  K-1) 

385 

700 

1978 

700 

903 

Table  3 

Specifications  of  the  LiMn204  battery  used  for  model  validation. 

Thickness  (Li) 

8.6  [mm] 

Width  (L2) 

115  [mm] 

Height  (I3) 

140  [mm] 

Weight 

316  [g] 

Nominal  voltage 

3.7  [V] 

Nominal  capacity 

8.5  [Ah] 

Specific  energy 

100  [Wh  kg-1] 

Energy  density 

230  [Wh  L"1] 

Maximum  cut-off  charge  voltage 

4.2  [V] 

Minimum  cut-off  discharge  voltage 

2  [V] 

Charge  working  temperature 

0-40  [°C] 

Discharge  working  temperature 

-20  to  55  [°C] 
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strongly  influenced  by  the  electrode  thickness  and  the  particle  size 
of  active  materials,  respectively.  All  of  these  affecting  factors  to  RH, 
direct  or  indirect,  will  be  discussed  in  details  in  the  following 
sections. 

3.1.  Model  validation 

To  validate  the  proposed  model,  a  discharge  test  is  performed  on 
a  commercially  available  Li-ion  battery  with  a  nominal  capacity  of 
8.5  Ah.  The  other  specifications  of  this  type  of  battery  are  given  in 
Table  3.  In  order  to  prevent  heat  transfer  to  the  surroundings,  a 
thermally  insulated  container  is  employed  in  the  test.  Both  the 
experimental  battery  and  the  simulated  battery  go  through  a  gal- 
vanostatic  discharge  process  at  a  1  C  rate  (i.e.,  8.5  A).  The  experi¬ 
ment  is  conducted  using  a  battery  analyzer  (ESI  PCBA  5010-4),  and 
the  temperatures  of  the  experimental  battery  are  collected  every 
60  s  by  a  K-type  thermocouple  attached  at  the  center  of  the  battery 
surface.  As  shown  in  Fig.  3,  the  simulation  result  is  well  consistent 
with  the  data  acquired  from  the  experiment,  which  demonstrates 
the  effectiveness  of  the  coupling  model  in  simulating  the  thermal 
response  associated  with  the  operation  of  Li-ion  battery. 

3.2.  Effect  of  discharge  rate 

To  investigate  the  influence  of  discharge  rate  on  the  total  heat 
generation  and  on  the  proportion  that  RH  accounts  for  the  total 
heat,  simulated  discharge  processes  at  rates  of  1  C  and  4  C  are 
carried  out,  through  which  the  RH  generation  rate  across  battery 
cell,  the  total  and  average  heat  generation  rates  with  and  without 
RH  included,  and  the  corresponding  temperature  responses  are 
inspected  as  follows. 

Fig.  4  shows  the  RH  generation  rate  across  the  cell  at  1  C 
discharge  rate,  in  which  the  current  collectors  are  not  involved 
since  no  RH  is  generated  on  these  two  foils.  According  to  Eq.  (13), 
RH  can  be  influenced  by  the  distribution  of  reaction  current  density 
as  well  as  the  variation  of  the  entropy  term  (dH/dT).  In  the  simu¬ 
lation,  the  cathode  is  thicker  than  the  anode  ( Fable  2),  which  makes 
the  distribution  of  reaction  current  density  is  relatively  uneven  on 
the  cathode.  Meanwhile,  during  the  course  of  discharge,  with  x  in 
Li*C6  decreases  from  0.56  to  0.04  and  y  in  LiyM^CH  increases  from 
0.15  to  0.7,  the  variation  of  the  entropy  term  given  in  Fig.  2 


Time  (s) 

Fig.  3.  Comparison  of  the  surface  center  temperatures  obtained  from  the  simulation 
and  experiment  at  1  C  discharge  rate  under  an  adiabatic  condition. 


Fig.  4.  Reversible  heat  generation  rate  across  the  cell  at  different  depths  of  1  C 
discharge. 


fluctuates  more  on  the  cathode  than  that  on  the  anode.  As  a 
result,  though  increases  remarkably  at  the  beginning,  the  RH  on  the 
anode  remains  stable  at  about  27,000  W  m“3  during  the  rest  period 
of  discharge,  while  the  variation  of  RH  generation  rate  at  the 
cathode  appears  more  complex  (between  -5000  and 


Fig.  5.  Total  heat  generation  rates  across  the  cell  at  different  depths  of  1  C  discharge  a) 
with  RH  source  considered;  and  b)  without  RH  source  considered. 
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5000  W  m~3).  And  in  general,  the  contribution  made  by  RH  to  the 
exothermic  behavior  of  the  battery  is  mainly  from  the  RH  generated 
on  the  anode. 

The  total  heat  generation  rate,  which  differs  dramatically  for  wr- 
model  and  wor-model,  is  also  presented  across  the  cell  at  various 
depths  of  discharge  (DOD)  to  identify  the  contribution  that  RH 
made  to  the  total  heat.  As  can  be  seen  from  Fig.  5,  by  including  RH, 
the  total  heat  generation  rate  at  the  anode  is  significantly  increased 
due  to  the  positive  RH  in  this  area,  with  the  average  values  more 
than  twice  that  of  the  model  neglecting  RH  at  most  of  the  selected 
DODs.  Whereas,  the  total  heat  generation  at  the  cathode  exhibits  a 
slight  decrease  after  introducing  RH  for  nearly  all  the  chosen  time 
points  except  1800  s  and  2400  s.  The  uneven  heat  distribution 
across  the  cell  should  also  be  noticed:  more  heat  is  generated  on  the 
anode  and  at  the  interfaces  between  electrode  and  separator,  which 
has  the  potential  to  cause  the  active  materials  close  to  separator  to 
degrade  fast  and  eventually  decreases  the  reliability  of  the  battery. 

Then  the  average  heat  generation  rates  and  the  temperature 
curves  of  the  3D  battery  models  with  and  without  RH  included  are 
depicted  in  Figs.  6a  and  7a,  respectively.  In  Fig.  6a,  the  solid  line 
indicates  the  heat  generation  including  RH,  whereas  the  dashed 
line  denotes  the  heat  generation  excluding  RH,  and  accordingly  the 
difference  between  these  two  lines  indicates  the  heat  generated 
from  the  RH  source.  When  ignoring  RH,  the  average  heat  genera¬ 
tion  rate  soars  in  the  first  few  seconds  and  remains  relatively  stable 
at  approximately  6000  W  m~3  for  most  of  the  time  prior  to  another 
increase  at  the  end  of  discharge,  which  is  believed  to  be  caused  by 
the  decrease  in  the  electrolyte  salt  conductivity  that  makes  the 
transportation  of  ions  more  difficult.  In  the  RH  inclusive  model, 
because  of  the  initially  negative  RH  generation  rate,  the  average 
heat  generation  rate  experiences  a  moderate  increase  initially  and 
then  a  continuous  increase  until  2500  s  when  the  value  of  the  solid 
curve  is  more  than  twice  that  of  the  dashed  curve.  It  eventually 
drops  as  a  result  of  the  decrease  in  the  RH  and  the  total  heat  on  the 


cathodes.  Correspondingly,  in  Fig.  7a,  the  RH  inclusive  model  pro¬ 
vides  a  lower  temperature  than  the  RH  exclusive  model  in  the 
beginning  of  discharge  (i.e.,  0-800  s),  but  a  faster  temperature 
increase  over  the  rest  period  of  discharge.  The  temperatures  of  the 
battery  models  with  and  without  RH  eventually  reach  322.5  K  and 
311.1  K,  respectively.  The  temperature  elevation  of  the  wr-model 
(24.35  K)  almost  doubles  that  of  the  wor-model  (12.95  I<). 

As  the  discharge  rate  increases  to  4  C,  the  RH  generation  also 
increases  in  view  of  Eq.  (13).  However,  the  proportion  of  RH  to  the 
total  heat  reduces,  as  indicated  by  the  decreased  relative  difference 
between  the  two  lines  presented  in  Fig.  6b  compared  to  Fig.  6a.  This 
is  mainly  due  to  the  dramatic  increase  of  the  ohmic  heat  which 
follows  an  i2R  type  of  heating  scheme.  As  a  result  of  the  reduction  of 
the  influence  of  RH  on  the  overall  heat  generation  across  the  cell, 
the  temperature  elevation  caused  by  RH  also  becomes  less  signifi¬ 
cant.  In  Fig.  7a,  at  the  end  of  discharge,  though  the  temperatures 
increase  to  as  high  as  340.8  K  and  331.6  K  for  models  with  and 
without  RH,  respectively,  the  ratio  of  the  temperature  changes 
(ATW/ATW0),  which  can  reveal  the  contribution  of  RH  source  made 
to  the  overall  temperature  rise,  is  calculated  to  be  1.28,  far  smaller 
than  the  1.88  for  the  1  C  discharge  case. 

3.3.  Effect  of  electrode  thickness 

As  a  main  geometrical  parameter  of  a  battery  cell,  the  thickness 
of  electrode  can  make  a  big  difference  in  the  electrochemical  and 
thermal  behaviors  of  Li-ion  battery  22].  Particularly,  a  variation  of 
the  electrode  thickness  will  cause  the  change  of  the  diffusion  dis¬ 
tance  inside  the  cell,  which  can  considerably  affect  the  amount  of 
ohmic  heat  produced,  and  so  do  the  percentages  of  different  types 
of  heat  that  constitute  the  overall  heat.  In  this  part,  the  thickness  of 
the  electrochemical  reaction  unit  (including  a  positive  electrode,  a 
negative  electrode  and  their  separator)  is  decreased  to  1  /4  of  the 
original  value  while  keeping  the  other  parameters  constant. 


Fig.  6.  The  average  heat  generation  rates  of  battery  models  with  (solid  line)  and  without  (dashed  line)  RH  included  for  the  a)  original  cell  at  1  C  discharge  rate,  b)  original  cell  at  4  C 
discharge  rate,  c)  cell  with  thin  electrode  at  1  C  discharge  rate,  d)  cell  with  small  particle  size  at  1  C  discharge  rate. 
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Fig.  7.  Temperature  curves  of  the  models  with/without  RH  included  for  a)  original  cell 
discharged  at  1  C  and  4  C  rates;  and  b)  adjusted  cells  discharged  at  1  C  rate. 


Considering  the  current  applied  on  the  current  collector  is  pro¬ 
portional  to  the  thickness  of  the  electrode,  it  decreases  as  the  same 
rate  as  the  electrode  thickness,  i.e.,  from  17.5  A  m-2  to  3.75  A  nrT2 
for  the  1  C  discharge  rate.  In  order  to  obtain  the  same  capacity  as 
the  original  battery,  3  times  more  reaction  cells  are  needed.  Overall, 
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Fig.  8.  Reversible  heat  generation  rate  across  the  cell  with  thin  electrodes  at  different 
depths  of  1  C  discharge. 
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the  thickness  of  the  newly  built  3D  thin-electrode  battery  thermal 
model  is  1.15  times  of  that  of  the  original  model. 

Fig.  8  shows  the  RH  generation  rate  across  the  adjusted  cell  at 
assigned  DODs.  Although  the  applied  current  decreases,  the  reac¬ 
tion  current  density  in  each  cell  remains  unchanged  owing  to  the 
decreased  electrode  thickness,  as  evidenced  by  the  similar  RH 
generation  rate  across  the  cell  before  and  after  changing  the  elec¬ 
trode  thickness. 

The  total  heat  generation  rates  of  the  RH  inclusive  and  exclusive 
models  are  given  in  Fig.  9a  and  b,  respectively.  As  the  electrode 
thickness  decreases,  the  total  heat  distribution  across  the  cell  ex¬ 
hibits  an  improved  uniformity,  which  can  benefit  the  long-term 
performance  and  stability  of  Li-ion  battery.  The  explanation  of 
the  improved  uniformity  is  that  during  the  operation  of  the  Li-ion 
battery,  the  movement  of  lithium  ions,  driven  by  the  potential 
and  concentration  gradients,  both  in  the  electrolyte  and  in  the 
active  particle  materials,  is  the  limited  step  (slower  than  the 
transportation  of  electrons).  Therefore  the  decrease  in  the  diffusion 
distance  can  facilitate  the  electrochemical  reactions  and  hence 
produces  a  better  distributed  current  and  reaction  depth,  a  reduced 
concentration  polarization,  a  less  ohmic  heat  generation  and  finally 
promotes  a  more  even  heat  distribution  across  the  battery  cell. 
Furthermore,  due  to  the  decrease  of  the  ohmic  heat,  the  proportion 
of  RH  to  the  overall  heat  increases,  as  is  evident  from  the  com¬ 
parison  made  between  Figs.  9  and  5. 
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Fig.  9.  Total  heat  generation  rates  across  the  cell  with  thin  electrodes  at  different 
depths  of  1  C  discharge  a)  with  RH  source  considered;  and  b)  without  RH  source 
considered. 
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Different  from  previous  models,  the  average  heat  generation 
rate  of  the  thin-electrode  battery  wr-model  (see  Fig.  6c)  experi¬ 
ences  a  sudden  drop  at  the  very  beginning  of  discharge,  which  is 
the  combined  effect  of  the  decreased  ohmic  heat  and  the  negative 
reversible  heat,  and  at  2500  s,  the  heat  generation  rate  of  the  wr- 
model  increases  to  almost  4  times  that  of  the  wor-model.  From 
the  big  difference  of  the  heat  generation  rates  of  the  wr-mode  and 
the  wor-model,  the  significance  of  RH  to  the  heat  generation  of  the 
battery  with  thin  electrodes  can  be  recognized. 

In  Fig.  7b,  an  alleviated  temperature  elevation  can  be  observed 
in  the  battery  with  adjusted  electrode  thickness.  Compared  to  the 
original  battery,  the  final  temperature  of  the  wr-model  decreases 
from  322.5  K  to  311  K,  which  not  only  reduces  the  security  risks  but 
also  enhances  the  well-being  and  cycle  life  of  the  battery.  Mean¬ 
while,  the  ratio  of  the  temperature  elevations  of  the  wr-model  and 
wor-model  becomes  significantly  bigger  in  the  adjusted  battery, 
2.45,  than  in  the  original  battery,  1.88.  In  light  of  the  large  ratio,  it  is 
fair  to  conclude  that  the  temperature  response  of  the  battery  will  be 
substantially  undervalued  by  excluding  RH  from  the  simulation, 
especially  for  the  battery  with  thin  electrochemical  reaction  cells. 

3.4.  Effect  of  active  material  particle  size 

The  solid  particle  size  of  electrode  active  materials  has  also  been 
acknowledged  as  an  important  factor  affecting  the  electrochemical 
processes  inside  the  Li-ion  battery  23,24].  The  decrease  in  particle 
size  can  facilitate  and  accelerate  the  processes  of  de-intercalation 
and  intercalation,  and  therefore  changes  the  time  and  depth  of 
the  electrochemical  reactions  taking  place  inside  the  battery  during 
discharge.  In  order  to  investigate  the  influence  of  particle  size  on 
the  thermal  behavior  of  the  battery,  in  this  part,  active  materials 
with  1  / 4  of  the  original  particle  size  for  both  anode  and  cathode  are 
used,  with  all  the  other  parameters  remain  the  same.  Also,  the 
adjusted  cell  with  the  new  particle  size  undergoes  a  1  C  discharge 
process. 

Prior  to  the  analysis  of  thermal  behavior,  a  comparison  made 
between  the  discharge  curves  of  the  battery  models  before  and 
after  modifying  particle  size  is  conducted.  As  depicted  in  Fig.  10,  the 
discharge  process  lasts  longer  in  the  battery  with  smaller  particles, 
which  allows  a  better  utilization  of  active  materials  and  hence 
provides  a  higher  capacity  output,  18.7  Ah  m~2,  as  compared  to  the 
17.5  Ah  m-2  of  the  original  battery.  The  explanation  is  based  on  two 


Fig.  10.  Discharge  curves  of  the  battery  models  with  original  particle  size  (dashed  line) 
and  1/4  of  the  original  particle  size  (solid  line). 


aspects:  first,  in  battery  with  smaller  particle  size,  the  difference 
between  the  Li-ion  concentration  at  the  surface  and  in  the  center  of 
the  electrode  particle  is  smaller  due  to  the  shorter  diffusion  dis¬ 
tances  for  both  de-intercalation  and  intercalation  reactions;  sec¬ 
ond,  during  discharge,  the  Li-ions  diffuse  from  the  center  of  the 
anode  particles  to  the  surface  to  react,  and  once  the  surface  Li-ion 
concentration  drops  to  zero,  the  reactions  stop.  Combining  these 
two  aspects,  at  any  time  point,  the  Li-ion  concentration  at  the 
surface  of  the  original  anode  particles  is  lower  than  that  of  the  small 
particles  due  to  the  longer  diffusion  distance.  And  eventually,  due 
to  the  lack  of  Li-ion  supply  on  the  anode  surface,  the  electro¬ 
chemical  reactions  in  the  original  battery  suffer  an  early  stop. 

Fig.  11  outlines  the  RH  generation  rate  across  the  reaction  cell 
with  modified  particle  size  at  assigned  DODs.  Comparing  with  Fig.  4 
(RH  of  the  original  cell),  the  decrease  in  particle  size  has  almost  no 
direct  impact  on  the  amount  of  the  RH  generated  except  promoting 
a  faster  response  and  a  better  material  utilization,  as  is  suggested  by 
the  intensive  endothermic  behavior  at  the  interfaces  between 
electrodes  and  separator  at  0  s.  And  at  any  time  point,  the  difference 
in  diffusion  distance  results  in  the  x  in  LixC6  lower  in  the  original 
particle  than  in  the  adjusted  particle,  and  thereby,  as  indicated  in 
Fig.  2,  the  RH  of  original  cell  always  corresponds  to  a  time  point  that 
is  ahead  of  the  time  point  to  which  the  RH  of  new  cell  corresponds. 
At  3600  s,  the  RH  across  the  anode  maintains  stable  at  approxi¬ 
mately  28,000  W  m”3  in  the  newly  built  cell,  whereas  it  distributes 
unevenly  in  the  original  cell  with  a  variation  of  10,000  W  nrT3  (from 
24,000  to  34,000  W  m“3),  which  is  attributed  to  the  non-uniform 
reaction  current  distribution  on  the  anode  as  well  as  the  uneven 
utilization  of  active  materials. 

The  active  material  particle  size  is  related  directly  to  the  gen¬ 
eration  of  reaction  heat,  which  is  less  produced  from  the  reaction 
over-potential  since  with  finer  particle  size  the  de-intercalation/ 
intercalation  reactions  in  the  battery  become  more  fluent.  Fig.  12 
presents  the  reaction  heat  generation  of  the  models  with  original 
and  decreased  particle  sizes.  As  shown,  the  reaction  heat  of  the 
anode  stays  at  around  7000  W  m~3  during  the  stable  stage  of 
discharge  in  the  original  cell,  while  it  reduces  to  only  about 
1700  W  m  3  for  the  cell  with  smaller  particle  size,  which  is  almost 
3/4  less,  and  the  change  is  basically  proportional  to  the  particle  size. 

As  a  consequence,  similar  to  the  case  of  thin-electrode  model, 
the  average  heat  in  the  model  with  small  active  material  particle 
has  a  higher  dependency  on  RH  than  in  the  original  one,  with  the 
proportion  of  RH  to  total  heat  at  2500  s  of  67.6%.  Moreover,  as  is 


Fig.  11.  Reversible  heat  generation  rate  across  the  cell  with  small  particles  at  different 
depths  of  1  C  discharge. 
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x  (m) 

Fig.  12.  Reaction  heat  generation  rate  across  the  cell  at  different  DODs  when  a)  1/4  of 
the  original  particle  size  and  b)  original  particle  size  are  used. 

represented  in  Fig.  7b,  the  temperature  curves  for  the  small-particle 
model  have  similar  trends  as  those  for  the  thin-electrode  model, 
with  the  final  temperatures  grows  higher  to  317.5  I<  (AT  =  19.35  K) 
and  306.5  K  (AT  =  8.35  K)  for  wr-model  and  wor-model,  respec¬ 
tively.  The  big  ratio  of  the  temperature  elevations  once  again  proves 
the  significance  of  taking  the  RH  into  consideration,  this  time  for 
cells  using  small  particle  size. 

To  give  a  brief  summary,  all  the  adjusted  parameters  and  the 
associated  results  are  listed  in  Table  4.  The  decreased  final  tem¬ 
peratures  for  the  thin-electrode  battery  and  for  the  small-particle 
battery  illustrate  an  improvement  in  the  thermal  performance  of 
these  batteries  with  respect  to  the  original  cell  under  1  C  discharge 
rate.  Note  that  the  high  temperature  of  the  original  battery  dis¬ 
charged  at  4  C  rate  should  be  properly  managed,  either  in  a  passive 
way  [25]  or  in  an  active  way  [26  ,  otherwise  it  can  degrade  the 
performance  and  reliability  of  Li-ion  battery,  or  even  bring  about 


Table  4 

Summary  of  the  parameter-varying  models  and  the  results  related  to  the  effect  of 
reversible  heat. 


Model 

C-rate 

Wo 

R/Rs 

qrev%  at  DOD  =  0.7 

ATw/ATwo 

Tfinal  (I<) 

1 

1 

1 

1 

57.6 

1.88 

322.5 

2 

4 

1 

1 

24.7 

1.28 

340.8 

3 

1 

0.25 

1 

72.6 

2.45 

311 

4 

1 

1 

0.25 

67.6 

2.32 

317.5 

Fig.  13.  Simulated  temperature  curve  of  the  original  battery  at  4  C  rate  with  forced  air 
convective  cooling. 

safety  issues,  such  as  thermal  runaway.  One  of  the  most  commonly 
used  active  cooling  approaches,  forced  air  cooling  method,  has  been 
adopted  in  the  simulation  to  control  the  temperature  rise  of  the 
original  battery  discharged  at  a  4  C  rate,  with  a  convective  coeffi¬ 
cient  value  of  40  W  m~2  K-1.  The  simulation  results  are  shown  in 
Fig.  13.  As  can  be  seen,  the  maximum  temperature  decreases  from 
340.8  I<  to  310.2  I<  when  the  forced  air  convective  cooling  is  applied. 

Fig.  14  indicates  the  combined  effects  of  the  C-rate  (1,  2,  4,  6,  8, 
10  C)  and  the  electrode  thickness  (0.1,  0.2,  0.4,  0.6,  0.8,  1  of  the 
original  cell)  on  the  results  of  Arw/ATwo  (due  to  the  similar  effects  of 
the  electrode  thickness  and  the  particle  size  on  thermal  response, 
the  figure  associated  with  the  impact  from  particle  size  is  not 
provided  here).  The  projections  of  the  obtained  data  points  are 
depicted  on  the  two  vertical  surfaces,  which  give  the  variation 
trend  of  the  ratio  of  the  temperature  rise  with  respect  to  the  C-rate 
and  the  relative  thickness.  The  big  ratios  shown  in  the  figure  reveal 
that  the  battery  with  thin  electrodes  and/or  discharging  at  low  C- 


Fig.  14.  Combined  effects  of  electrode  thickness  and  discharge  rate  on  the  ratio  of 
temperature  elevations  in  models  including  and  excluding  reversible  heat. 
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rates  can  produce  significant  proportion  of  RH  to  the  total  heat,  and 
although  this  proportion  decreases  as  the  electrodes  become 
thicker  and/or  the  discharge  rate  becomes  higher,  the  contribution 
of  RH  should  be  taken  into  account  for  discharge  processes  at  10  C 
or  lower  rates.  Furthermore,  this  diagram  is  expected  to  be  used  for 
correcting  the  final  temperatures  obtained  in  the  previous  and 
future  thermal  analyses  of  Li-ion  batteries  that  do  not  include  the 
reversible  heat. 

4.  Conclusion 

In  this  paper,  a  reversible-heat-included  Li-ion  battery  model 
consisting  of  a  ID  electrochemical  part  and  a  3D  thermal  part  is 
developed  and  used  to  simulate  several  galvanostatic  discharge 
processes  on  batteries  with  different  physical  properties.  By 
comparing  the  results  obtained  from  the  wr-model  and  from  the 
wor-model,  a  better  understanding  of  the  importance  of  RH  to  the 
overall  thermal  performance  of  Li-ion  battery  under  different 
conditions  is  achieved,  based  on  which  a  correction  diagram  is 
proposed  to  modify  the  results  for  the  thermal  analyses  on  Li-ion 
batteries  that  do  not  include  RH  at  first. 

RH  is  non-ignorable  in  battery  discharged  at  low  rates,  but  its 
contribution  to  the  total  heat  decreases  as  a  result  of  the  dramatic 
increase  of  the  ohmic  heat  at  high  rates.  Battery  with  thinner 
electrodes  relies  heavier  on  the  RH  than  that  with  thicker  elec¬ 
trodes.  Instead  of  the  change  of  the  RH  itself  with  the  reduction  of 
the  thicknesses  of  electrodes,  the  considerable  decrease  in  ohmic 
heat  generation  takes  the  main  responsibility,  which  indirectly 
increases  the  contribution  of  RH  to  the  total  heat  generation. 
Similarly,  the  effect  of  RH  is  greater  in  battery  with  smaller  active 
material  particle  size,  as  a  result  of  the  decreased  reaction  heat. 

Moreover,  a  strong  correlation  between  the  battery  temperature 
and  the  battery  cell  characteristics  is  demonstrated.  It  is  shown  that 
the  decrease  in  any  of  these  two  parameters,  i.e.,  electrode  thick¬ 
ness  and  active  material  particle  size,  can  reduce  the  battery  tem¬ 
perature  during  operation  and  thereby  improving  the  performance, 
reliability  and  even  safety  of  Li-ion  battery. 

Overall,  the  significance  analysis  of  reversible  heat  is  crucial 
when  precise  information  about  the  heat  generation  and  the  heat 
distribution  across  a  battery  cell  is  demanded  and  should  be  con¬ 
ducted  case-by-case  since  the  proportion  of  RH  varies  dramatically 
with  the  applied  currents  and  with  the  characteristics  of  the  elec¬ 
trochemical  battery  under  study. 
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F  Faraday's  constant,  96,487  C  mol-1 

h  convective  heat  transfer  coefficient,  W  nrT2  K_1 

ie  solution  phase  current  density,  A  cm-2 

is  solid  phase  current  density,  A  cm-2 

jii  reaction  current  density,  A  cm-3 

k  conductivity  of  different  components  inside  cell, 

W  m  1  K"1 

/<eff  effective  electrolyte  conductivity,  S  m_1 

lc  thickness  of  copper  foil,  cm 

ln  thickness  of  negative  electrode,  cm 

ls  thickness  of  separator,  cm 

lp  thickness  of  positive  electrode,  cm 

la  thickness  of  aluminum  foil,  cm 

lv  thickness  of  components  in  experimental  cell,  cm 

l  thickness  of  entire  electrochemical  cell,  cm 

L\  thickness  of  experimental  battery,  cm 

L2  width  of  experimental  battery,  cm 

L3  length  of  experimental  battery,  cm 

n  number  of  electrons  transferred  in  the  electrode  reaction 

q  heat  generation  rate,  W  m-3 

qr  reaction  heat  generation  rate,  W  m-3 

qrev  reversible  heat  generation  rate,  W  m~3 

qj  ohmic  heat  generation  rate,  W  m~3 

q  charge  density  on  the  double  layer  surface,  C  m-2 

r  radial  coordinate  along  active  material  particle,  cm 

R  universal  gas  constant,  8.3143  J  mol-1  K-1 

Si  stoichiometric  coefficient 

t  time,  s 

t°  transference  number 

T  temperature,  K 

U  open  circuit  potential,  V 

dU/dT  entropy  change  term,  mV  K-1 

Dref  open-circuit  potential  of  an  electrode  reaction,  V 

x  x  coordinate 

x  negative  electrode  solid  phase  stoichiometry 

y  y  coordinate 

y  positive  electrode  solid  phase  stoichiometry 

z  z  coordinate 


Greek  symbols 

e  volume  fraction  or  porosity  of  a  phase 

ri  surface  over-potential  of  an  electrode  reaction,  V 

k  reaction  rate  constant 

a  conductivity  of  solid  active  materials  in  an  electrode, 

S  cirT1 

p  density  of  the  battery,  kg  m~3 

v+  number  of  cations  of  anions  into  which  a  mole  of 

electrolyte  dissociates 
$  electrical  potential  in  a  phase,  V 


Nomenclature 

as  active  surface  area  per  electrode  unit  volume,  cm2  cm-3 

c  volume-averaged  concentration  of  lithium  in  a  phase, 

mol  cm-3 

cs,max  maximum  concentration  of  lithium  in  solid  phase, 
mol  cm-3 

cs>o  initial  concentration  of  lithium  in  solid  phase,  mol  cm-3 
Cdi  double  layer  capacitance,  F  m~2 

Cp  specific  heat  capacity  of  component  inside  cell,  J  kg-1 1<-1 

D  diffusion  coefficient  of  lithium  species,  cm2  s 

Ds  diffusion  coefficient  of  lithium  in  the  electrode  particles, 

cm2  s 

e~  electrons 


Subscripts 
amb  ambient 

c  convection 

e  electrolyte  phase 

f  conductive  filler  additive 

i  ith  component 

p  polymer  phase 

ref  with  respect  to  a  reference  state 

s  solid  phase 


Superscripts 

eff  effective  value 


432 


R.  Zhao  et  al.  /  Journal  of  Power  Sources  266  (2014)  422—432 


References 

[1  ]  H.  Li,  X.  Huang,  L.  Chen,  Z.  Wu,  Y.  Liang,  Electrochem.  Solid-State  Lett.  2  (1999) 
547-549. 

[2]  C.C.  Lin,  H.C.  Wu,  J.P.  Pan,  C.Y.  Su,  T.H.  Wang,  H.S.  Sheu,  N.L.  Wu,  Electrochim. 
Acta  101  (2013)  11-17. 

[3]  L.  Xia,  D.  Wang,  H.  Yang,  Y.  Cao,  X.  Ai,  Electrochem.  Commun.  25  (2012) 
98-100. 

[4]  H.  Kawai,  M.  Nagata,  H.  Tukamoto,  A.N.  West,  J.  Power  Sources  81—82  (1999) 
67-72. 

[5]  T.M.  Bandhauer,  S.  Garimella,  T.F.  Fuller,  J.  Electrochem.  Soc.  158  (2011) 
R1-R25. 

[6]  W.B.  Gu,  C.Y.  Wang,  ECS  Proc.  99  (25)  (2000)  748-762. 

[7]  B.  Wu,  V.  Yufit,  M.  Marinescu,  G.J.  Offer,  R.F.  Martinez-Botas,  N.P.  Brandon, 
J.  Power  Sources  243  (2013)  544—554. 

[8]  V.  Srinivasan,  C.Y.  Wang,  J.  Electrochem.  Soc.  150  (2003)  A98— A106. 

[9]  Y.  Reynier,  J.  Graetz,  T.  Swan-Wood,  P.  Rez,  R.  Yazami,  B.  Fultz,  Phys.  Rev.  B  70 
(2004)  174304. 

[10]  Y.  Reynier,  R.  Yazami,  B.  Fultz,  J.  Electrochem.  Soc.  151  (2004)  A422— A426. 

[11]  S.W.  Kim,  S.I.  Pyun,  Electrochim.  Acta  46  (2001)  987—997. 

[12]  R.  Williford,  V.  Viswanathan,  J.G.  Zhang,  J.  Power  Sources  189(2009)  101—107. 


[13]  D.H.  Jeon,  S.M.  Baek,  Energy  Convers.  Manage.  52  (2011)  2973—2981. 

[14]  K.  Smith,  C.Y.  Wang,  J.  Power  Sources  160  (2006)  662—673. 

[15]  G.  Sikha,  R.E.  White,  B.N.  Popov,  J.  Electrochem.  Soc.  152  (2005) 
A1682— A1693. 

[16]  M.  Doyle,  J.  Newman,  A.S.  Gozdz,  C.N.  Schmutz,  J.M.  Tarascon,  J.  Electrochem. 
Soc.  143  (1996)  1890-1903. 

[17]  K.E.  Thomas,  C.  Boguta,  J.  Newman,  J.  Electrochem.  Soc.  148  (2001) 
A570-A575. 

[18]  S.  Al  Hallaj,  R.  Venkatachalapathy,  J.  Prakash,  J.R.  Selman,  J.  Electrochem.  Soc. 
147  (2000)  2432-2436. 

[19]  P.  Taheri,  M.  Yazdanpour,  M.  Bahrami,  J.  Power  Sources  243  (2013)  280—289. 

[20]  S.C.  Chen,  C.C.  Wan,  Y.Y.  Wang,  J.  Power  Sources  140  (2005)  111—124. 

[21]  U.S.  Kim,  C.B.  Shin,  C.S.  Kim,  J.  Power  Sources  189  (2009)  841—846. 

[22]  H.  Zheng,  J.  Li,  X.  Song,  G.  Liu,  V.S.  Battaglia,  Electrochim.  Acta  71  (2012) 
258-265. 

[23]  C.H.  Lu,  S.W.  Lin,  J.  Power  Sources  97-98  (2001)  458-460. 

[24]  Y.  Chen,  K.  Xie,  Y.  Pan,  C.  Zheng,  J.  Power  Sources  196  (2011)  6493—6497. 

[25]  R.  Zhao,  S.  Zhang,  J.  Gu,  J.  Liu,  S.  Carkner,  E.  Lanoue,  J.  Power  Sources  255 
(2014)  29-36. 

[26]  H.  Park,  J.  Power  Sources  239  (2013)  30—36. 


